Two-photon exchange in elastic electron-proton scattering : QCD factorization 

approach 
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We estimate the two-plioton excliange contribution to elastic electron-proton scattering at large 
momentum transfer . It is shown that the leading two-photon exchange amplitude behaves as 
l/Q* relative to the one-photon amplitude, and can be expressed in a model independent way in 
terms of the leading twist nucleon distribution amplitudes. Using several models for the nucleon 
distribution amplitudes, we provide estimates for existing data and for ongoing experiments. 
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Elastic electron-nucleon scattering in the one-photon 
(I7) exchange approximation is a time- honored tool for 
accessing information on the structure of the nucleon. 
Precision measurements of the proton electric to mag- 
netic form factor ratio at larger using polarization 
experiments P, 0, Q have revealed significant discrepan- 
cies in recent years with unpolarized experiments using 
the Rosenbluth technique y. As no experimental flaw 
in either technique has been found, two-photon (27) ex- 
change processes are the most likely culprit to explain 
this difference. Their study has received a lot of attention 
lately, see @ for a recent review (and references therein), 
and [6[ for a recent global analysis of elastic electron- 
proton (ep) scattering including 2j corrections. In this 
work we calculate the leading in behavior of the elas- 
tic ep scattering amplitude with hard 27 exchange. 

To describe the elastic ep scattering, l{k) + N{p) — > 
l{k') + N{p'), we adopt the definitions : P = {p + p')/2, 
K ^ {k + k')/2, q = k-k' = p' -p, and choose = -q^ 
and V — K ■ P as the independent kinematical invariants. 
Neglecting the electron mass, it was shown in [7] that 
the T-matrix for elastic ep scattering can be expressed 
through 3 independent Lorentz structures as : 
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where e is the proton charge, M is the proton mass, 
h = ±1/2 is the electron helicity and Xn (A'^) are 
the helicities of the incoming (outgoing) proton. In 
Eq. ^ , Gm , F2, F3 are complex functions of v and . 
To separate the I7 and 27 exchange contributions, it 
is furthermore useful to introduce the decompositions : 
Gm = Gm + SGm, and ^ F2 + SF2, where GAf(F2) 
are the proton magnetic (Pauli) form factors (FFs) re- 
spectively, defined from the matrix element of the elec- 
tromagnetic current, with Gm(0) = ^ip — 2.79 the pro- 
ton magnetic moment. The amplitudes F3, SGm and dF2 
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FIG. 1: Typical graph for the elastic ep scattering with two 
hard photon exchanges. The crosses indicate the other possi- 
bilities to attach the gluon. The third quark is conventionally 
chosen as the d— quark. There are other diagrams where the 
one photon is connected with u— and d— quarks. We do not 
show these graphs for simplicity. 



originate from processes involving at least 27 exchange 
and are of order (relative to the factor in Eq. ([T|)). 

The leading perturbative QCD (pQCD) contribution 
to the 27 exchange correction to the elastic ep ampli- 
tude is given by a convolution integral of the proton 
distribution amplitudes (DAs) with the hard coefficient 
function as shown in Fig. [TJ In the hard regime, where 
Q^,s ^ M^, we calculate the amplitude in the Breit 
system, where the initial and final proton momenta cor- 
respond to two opposite light-like directions : 

TL 

p ~ Q-, with n = (1,0,0,1), 
p' ^ g-, with n= (1,0,0,-1), 

with {n • n) — 2. The lepton kinematics are given by : 

71 — Ti —71 71 

k^CQ^- CQ 2 + k^, k' = -CQ- + + ^-L' 

where, at large Q^, C can be determined from s ~ CQ^i 
and u ~ — CQ^i with ^ = 1 — and C ^ 1- Further- 
more, the transverse vector in the lepton kinematics is 
determined from : fc^ = — CCQ^- 

Let us now consider the proton matrix element which 
appears in the graph of Fig. [T] Following the notation 
from the proton matrix element is described at lead- 
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ing twist level by three nucleon DAs as : 
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with light-cone momentum p+ = Q, where C is charge 
conjugation matrix : C~^^^C = — 7J, and where X = 
{A, V, T} stand for the nucleon DAs which are defined by 
the light-cone matrix element : 

Xia,,\p+) I d[x,] e-'^P^(^^"^')x(x,), (3) 



with 

d[xi] = dxidx2dx^5{l — Xj). 
In general, the following properties are valid 

V{xi,X2,X3,) = V{x2,Xi,Xz), 
A{xi,X2,X3,) = -A{x2,Xi,Xz), 

T{xi,X2,X3) = i[y-A](l,3,2) + i[y-A](2,3,l), 

i.e. we have only two independent functions. 

In the large limit, the pQCD calculation of Fig. [T] 
involves 24 diagrams, and leads to hard 27 corrections to 
SGm, and i^/APF^, which are found as : 

oGm — 
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X / d[yi] d[xi 



4x2 y2 
D 

X {Qu [{V + A'){V + A)+ AT'T] (3, 2, 1) 
+QuQd [{V + A')iV + A)+ AT'T] (1, 2, 3) 



+QuQd2[V'V + A]{l,i,2)} , 



(4) 



and 

V 



2(3:2 j/2 + ^2 2/2) 



X {Qu' [(1^' + A'){V + A)+ AT'T] (3, 2, 1) 
+QnQd [(V^' + A'){V + A)+ AT'T] (1, 2, 3) 
+QuQd2[V'V + A'A]{l,i,2)}, (5) 

with quark charges Qu ~ +2/3, Qd — —1/3, aem — 
e^/(47r), Usip?) is the strong couphng constant evaluated 
at scale , and the denominator factor D is defined as : 

D = {yiy2y2) {X1X2X2) [2:2 C + 2/2 C - 2^22/2 + ie] 

X [X2C + 2/2C - 2:2^2 + is] ■ (6) 



The unprimed (primed) quantities in Eqs. (|4l [5]) refer to 
the DAs in the initial (final) proton respectively. Eqs. ([H 
[S]) are the central result of the present work. One no- 
tices that at large Q^, the leading behavior for 6Gm and 
v/M'^F^ goes as l/Q'*. In contrast, the invariant 5F2 is 
suppressed in this limit and behaves as l/Q^. 

It is interesting to point out that the scaling behavior 
for the 27 amplitude obtained in the present calculation 
differs from the handbag calculation of [^, [l3] ■ Whereas 
the present calculation gives a model independent leading 
behavior of l/Q** for the 27 amplitude relative to the I7 
amplitude, the behavior of the 27 amplitude within 
the handbag calculation depends on the specific modeling 
of the generalized parton distributions. As an exam ple , 
the modified Regge parameterization considered in (9l.Tlo| 
leads to a calculation which is of higher twist compared 
to the leading pQCD calculation considered here. 

To evaluate the convolution integrals in Eqs. (HlISl), we 
need to insert a model for the nucleon twist-3 DAs, V , 
A, and T. The asymptotic behavior of the DAs and their 
first conformal moments were given in 



as 



V{x^) ~ l2QxiX2X:ifN[l + r+{l-ix^)] 
A{xi) ~ l2QxiX2X'ifN r^{x2 - xi), 



T{xi) ~ l2QxiX2X-ifN 



1 + - (r_ - r+) (1 - 3x3) 



(7) 



and depend on three parameters : /at, r_ and r+. In 
this work, we will provide calculations using two models 
for the DAs that were discussed in the literature. The 
corresponding parameters (at ^ = 1 GeV ) are given in 
Table |T1 and are compared with recent lattice QCD cal- 
culations (QCDSF flJl), extrapolated to the chiral limit. 
One notices that the parameters r_ and r+ in the BLW 
model for the proton DA are totally compatible with the 
lattice results, whereas the value of the overall normal- 
ization /jv for the lattice DA is about 2/3 smaller than 
the BLW value. Below, we will provide calculations using 
the models COZ and BLW, and note that a good esti- 
mate using the lattice DA can be directly obtained from 
our figures by scaling the BLW result by a factor « 2/3. 





(10-^ GeV^) 


r- 


r+ 


COZ fll] 


5.0 ±0.5 


4.0 ± 1.5 


1.1 ±0.3 


BLW [12] 


5.0 ±0.5 


1.37 


0.35 


QCDSF fl3] 


3.23 
±0.06 ± 0.09 


1.06 
±0.09 ±0.31 


0.33 
±0.03 ±0.11 



TABLE I: Parameters entering the proton DA (at ^ = 1 GeV) 
for two models (COZ, BLW) used in this work. For compari- 
son we also show a recent lattice evaluation (QCDSF). 

We next calculate the effect of hard 27 exchange, given 
through Eqs. ([H [5]), on the elastic ep scattering observ- 
ables. The general formulas for the observables including 
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the 27 corrections 5Gm, SF2, and v/M^F^ were derived 
in , to which we refer for the corresponding expressions. 
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FIG. 2: Rosenbluth plots for elastic ep scattering: an divided 
by /Xp/(1 + gVO.Tl)*. Dashed (blue) curves : I7 exchange, 
using the Gep/Gmp ratio from polarization data 0, 0: [3- 
Sohd red (dotted black) curves show the effect including hard 
27 exchange calculated with the BLW (COZ) model for the 
proton DAs. The vertical dotted line shows the boundary 
where the Ihs of Eq. ([8]) is 0.5. The data are from Ref. i_4]. 

In Fig. O we calculate the reduced cross section 
as a function of the photon polarization parameter e 
and different values of Q^. In the I7 exchange, aji = 
Gm{Q^) + e/r G£;(Q2), with r = Q^/iAM"^), and the 
Rosenbluth plot is linear in e, indicated by the dashed 
straight lines in Fig. [21 The effect including the hard 
27 exchange is shown for both the COZ and BLW mod- 



els of the proton DAs. One sees that including the 27 
exchange changes the slope of the Rosenbluth plot, and 
that sizeable non-linearities only occur for e close to 1. 
The inclusion of the hard 27 exchange is able to well de- 
scribe the dependence of the unpolarized data, when 
using the polarization data P, 0, S] for the proton FF ra- 
tio Gep/Gmp as input. Quantitatively, the COZ model 
for the nucleon DA leads to a correction about twice as 
large as when using the BLW model. The question arises 
as to the applicability of the hard scattering calculation 
for the values of the data shown in Fig. [21 On the one 
hand, the argument of the running coupling as is defined 
by the renormalization scale Q, which should be suf- 
ficiently large to validate a pQCD calculation. On the 
other hand, contains a logarithmic singularity when 
e — » (i.e. when ^ — > 1). Therefore, the applicability of 
the hard description is restriced by the condition : 

a,(Q2) |ln(l-l/C)| «L (8) 

In Fig [21 we indicate a boundary by dashed vertical lines, 
where q;5(Q^) |ln(l — 1/C)| < 0.5, which corresponds 
with £ > 0.25. We hke to note here that in contrast to 
the pQCD treatment of the proton FFs, which requires 
two hard gluon exchanges, the 27 correction to elastic ep 
scattering only requires one hard gluon exchange. One 
therefore expects the pQCD calculation to set in for 
values in the few GeV^ range, which is well confirmed by 
the results shown in Fig. [21 

The real part of the 27 exchange amplitude can be ac- 
cessed directly as the deviation from unity of the ratio 
of e^/e~ elastic scattering. The precision of past ex- 
periments performed at SLAC [3], was not sufficient to 
see a clear deviation from unity over a large range in 
e. Presently, several new experiments are planned or are 
underway at VEPP-3 JLab/CLAS 0, and Olym- 
pus@DESY to make precision measurements of the 
e^/e" ratio in elastic scattering off a proton. In Fig. [3l 
we show the predictions for the (Te+p/cye-p ratio for differ- 
ent values ofQ^ and e planned by the 01ympus@DESY 
experiment [17||. In order to make a comparison with our 
pQCD calculations, we only show kinematics for which 
> 2 GeV^. We also made an estimate of the theoret- 
ical uncertainty of our predictions. For this purpose, we 
varied the normalization scale /i (entering as) and the 
normalization Jn (in units 10"'^ GeV^) for both COZ 
and BLW models of the proton DAs over the ranges : 
< fx^ < Q^, and 4.5 < Jn < 5.5. Furthermore, for 
the COZ model, we also varied the parameters r+ and r_ 
in the range : 0.8 < r+ < 1.4, and 2.5 < r_ < 5.5. The 
resulting theoretical error estimate is also indicated on 
Fig. [3l The 01ympus@DESY experiment aims at a sta- 
tistical precision of the e^/e" cross section ratio of better 
than one percent for an average = 2.2 GeV^. Our cal- 
culations predict a deviation from unity for these ratios 
in the range 2.5% (BLW) to 5% (COZ), allowing for an 
unambiguous test with the upcoming measurements. 
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FIG. 3: Predictions for the ratio cr^+p/o"e-p for different val- 
ues of (shown by numbers in GeV^) and e planned by the 
01ympus@DESY [l7|l experiment (we only show kinematics 
for which > 2 GeV^). The upper blue (lower red) points 
correspond with the COZ (BLW) models. The error-bars 
show the theoretical uncertainties due to the parameters of 
the proton DAs and the running coupling scale, as described 
in the text. For comparision, we also show the theoretical 
curves for the average values — 2.4 GeV'^ (dashed curves) 
and = 3.25 GeV^ (dotted curves) for both the COZ (two 
upper curves) and BLW (two lower curves) models. 
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FIG. 4; The ratio Ps/Pi as a function of e for = 2.5 GeV^ . 
The horizontal curve is the result of the I7 calculation. The 
dotted black (solid red) curves correspond to the I7 -I- 27 cal- 
culations using the COZ (BLW) models for the proton DAs. 
The data point is from the JLab/Hall A experiment P, [!]■ 



In Fig. m we show the ratio of the proton recoil 
polarization components Ps/Pi as a function of e for 
= 2.5 GeV^. The I7 prediction is given by the hor- 
izontal line. Our I7 -I- 27 prediction yields a negative 
correction which increases with decreasing e. Around 
s = 0.3 it reduces the PjPi ratio by 2 % for the BLW 
model, and by around 4 % for the COZ model. The 
JLab/Hall A experiment has measured this ratio at 
a large e value around 0.85. A new JLab/Hall C experi- 
ment jl8| . which is currently under analysis, has recently 
measured this ratio at ~ 2.5 GeV^ for three e values 
between 0.15 and 0.8. The expected experimental preci- 
sion of around 1% for this ratio, will therefore allow to 
test our predictions, which are in the 2 — 5 % range. 

The observables discussed above test the real parts 
of the 27 amplitudes. The imaginary parts of Eqs. ([H 



[5]), arising from lepton propagator singularities, can be 
tested by polarizing the target or recoiling proton per- 
pendicular to the scattering plane. 

In summary, we calculated the leading in behavior 
of the 27 exchange contribution to elastic ep scattering. 
It was found that the leading 27 amplitude is given by 
processes involving one hard gluon exchange, resulting 
in a behavior of the 27 amplitude relative to the 

I7 amplitude. We expressed the leading 27 amplitude 
in terms of the leading twist nucleon DAs. Using two 
models for the nucleon DAs, we found that, for in 
the few GeV^ range, these calculations can quantitatively 
explain the slope of the Rosenbluth plot when using the 
Gep/Gmp polarization data as input. Furthermore, we 
have shown that ongoing and planned clastic ep scatter- 
ing experiments both for the s dependence of the recoil 
polarization ratio Ps/Pi as well as for the e+/e~ ratio, 
have the precision to test our predictions. 
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and the U.S. DOE under contract DE-FG02-04ER41302. 
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